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(~| Abstract 

-(— > 

^ We consider 2-diniensional Schrodinger operator with the non- 

degenerating magnetic field and we discuss spectral asymptotics with 
,— I the remainder estimate o{fi^^ h^^) or better. 

K*" We also consider 3-dimensional Schrodinger operator with the 

QQ non-degenerating magnetic field and we discuss spectral asymptotics 

■^ with the remainder estimate o{iJ,~^ h~^) or better. 

^^ We also consider generalized Schrodinger-Pauli operator in the 

^O same framework albeit with fih > 1 and derive spectral asymptotics 

with the remainder estimate up to 0{h \ log h\) and with the princi- 
pal part X fih~'^. 
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^ Introduction 



Our goal is to derive spectral asymptotics of 2, 3-dimensional Schrodinger 
operator 

(0.1) A = J2Pjg^''Pk+V, with Pj = hDj - fiVj 

j.k 

near the boundary where 
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(0.2) g-''^ = g'^J , Vj and V are real-valued functions, h G (0, 1], /i G [1, oo). 

and in 6(0, 1) C X the following conditions are fulfilled: 

(0.3)i_3 |DV1 < c, \D'^Vj\<c, |D"\/|<c ya:\a\<K, 

(0.4) tQ<^gJ^r]jT]i,-\T]\-^ <c V?7GM^\0 VxG 6(0,1) 

and unless opposite is specified, we always assume that 

(0.5) X 3 6(0,1). 

So we basically want to generalize results of sections 13.3-13.5 of Chap- 
ter 13 [lvr2]-'^\ We assume that condition 

(0.6) F > eo 

where F is the scalar intensity of magnetic field defined by (13.1.13)23, as 
d = 2,3 respectively. 

What we want is to analyze how improved dynamics results in the sharper 
remainder estimates. 

Plan of the article 

We start with the two-dimensional case when singularities propagate along 
magnetic drift lines and then proceed to the three-dimensional case when 
singularities propagate along magnetic lines which could however have "side- 
drift". 



^' This article is a rather smaU part of the huge project to write a book and is just part 
of the section 13.6 consisting entirely of newly researched results. Chapter 13 corresponds 
to Chapter 6 of its predecessor V. Ivrii [Ivrl]. External references by default are to [Ivr2]. 
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In the main section 3 we consider supersharp asymptotics as c/ = 3 with 
the remainder estimate 0{fi~^h^^) as fih < 1 and 0{h~^) as ^h > 1. 

1 Case d = 2 

1.1 Classical dynamics and heuristics 

Let us start from the classical dynamics. We assume that potential is truly 
generic i.e. no special conditions with u are assumed. To make things 
simpler we assume at the moment that g-''^ = 5jk and F = const. 

Looking at the canonical form we conclude that there is fast circular 
movement in {fiX2,C,2) and there is a slow drift movement in {xi, fi~^C,i)- 
More precisely, consider first two terms of the canonical form: 

(LI) wo(xi,/^-'6) + wi(xi,/i-^ei)(A2 +6) + 0{fi~^) 

and we understand that along {fiX2, ^2) there is a circular movement with 
period 7"(/x) = irfi^^cu^^ + 0{fi^^) and along {xi,fi~^^i there is a drift with 
the velocity fi^^ H^^^j^^^^k with Hamiltonian calculated with respect to (xi,,^i) 
only and with k = [t — UJo)ui plugged after (as we are interested at energy 
level 0). 

Going back with symplectomorphism '^ we conclude that modulo 0{n^^) 

(1.2) 7(/i) = 7r/i-^F-^ 
and the drift speed is 

(1.3) v,,f, = /.-^(V(^))^ 

which defines magnetic drift lines which are level lines of (V — t)F^^ "^K 
Remark 1.1. Note that in an invariant form 

(L4) F = ^|Fi2|, Fi2 = 9iV2-92Vi, g = (gV'-grV)~' 



'' In the special case F — const magnetic drift lines are level lines of V . 
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and magnetic drift is described by 



(1.5) 

where 
(1.6) 

is a pseudo-scalar^\ 



dXi ■'■ a ^ ^^2 

~dt^7/^ ^7' ~dt 
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Remark 1.2. Magnetic drift lines are terminated in the critical points of 
(\/-r)F-i4). 

If we are interested in the time T ^ /i we need to follow magnetic drift 
line either close to termination point, or to large distance, or assume that 
it is closed smooth line. However we need to remember that Hamiltonian 
trajectories follow drift lines albeit do not coincide with them! 

Example 1.3. Let g^^ = bjk, F = 1. 

(i) As V = X2 magnetic drift lines are just straight lines X2 = const with no 
termination points and Hamiltonian trajectories are cycloids of figure 1(b). 





(a) Periodic 



(b) Drift 



Figure 1: Trajectories for 2-dimensional model operator: 
a) Unperturbed b) Perturbed by a constant electric field 



(ii) Consider V = ±(x^ + x|). Then magnetic drift lines are circles x^ + 
x| = const. On the other hand, Hamiltonian trajectories are cycloids with 

'^' I.e. it changes sign as we change orientation of the coordinate system. 
^^ In the special case F — const those are critical points of V. 
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(a) Helix (b) Deformed helix 

Figure 2: Trajectories for 3-dimensional model operator: 

a) Unperturbed b) Perturbed by a constant electric field parallel to magnetic 

lines: non-uniform movement along magnetic lines; drift is negligible in 

comparison 



frequencies 2/i(l + 0{fi ^) and 2/i ^ + (1 + 0{fi ^) and they are periodic 
only as /i = /i/( = /c2(l + 0{k~^) with k G Q+. See figure 3(a), (b). 

(iii) Consider V = xiX2. Then magnetic drift lines are hyperbolas xiX2 = 
const. On the other hand, Hamiltonian trajectories are shown on fig- 
ure 3(c),(d), one which terminates in x = and another which bypasses it. 



Our conjecture is that we can derive remainder estimate o(^fj,^^h^^ + l) 
in all the cases when drift lines are not closed; if drift lines are closed we 
can derive this remainder estimate for /i except belonging to rather thin 
exceptional set. 

One should notice that magnetic drift lines escaping to infinity are rather 
dangerous as then V does not grow which may indicate presence of the 
essential spectrum. 

2 Case d = 3 

Now we consider 3-dimensional case. 
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(a) Periodic 



(b) Non-periodic 





(c) Termination 



(d) Non-termination 



Figure 3: Trajectories for 2-dimensional model operator to example 1.3(ii) 
(in (a),(b)) and 1.3(iii) (in (c),(d)) 



2.1 Preliminaries and improved reduction 

It is well known from physics that the classical three-dimensional particles 
in a strong magnetic field move almost along magnetic lines. So it is not 
surprising that in our analysis magnetic lines are also important. Let us 
recall that a magnetic line of the magnetic field F = [F^, F^, F'^) is an 
integral curve x = x[t) of the field F~^f. Thus, x[t) is given by the system 

dxj 2 



(2.1) 



dt 



-.F^ a = 1.2, 3) 
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where 

(2.2) ^ = (E SjkPF') ' = 1{J2 S^'g'^'FjiF,^) ^ 

j,k j,k,l,m 



is the scalar intensity; let us introduce the corresponding flow ct>t on X. 
It follows from the definition of F that this field is solenoidal, i.e., 

(2.3) 5^9;(Fv^) = 



and therefore 

(2.4) Magnetic flow (t>t preserves the density ^Jgdx on X. 
Moreover, 

(2.5) Magnetic flow 0f can be lifted to the Hamiltonian flow ^i on A which 
is given by the Hamiltonian 

(2.6) ^E^'fe-'^^iW)- 

Remark 2.1. We could deflne the flow <t>t by the system 

(2.1)' $ = ?"^' a = 1.2, 3) 

with functions u and uj~-^ regular along the trajectory; then 4>f would be 
given by the Hamiltonian 



(2.6)' ^E'''(0--"^-(^)) 



J 

-1, 



and the invariant density would be to F^/gdx. 

Remark 2.2. Even from the heuristic point of view there is an essential 
difference: 

(i) Magnetic lines depend only on F while magnetic drift lines depend on V 
as well. In order to the get magnetic drift lines at the energy level r one 
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should replace V hj V — t. Then, for F ^ const the magnetic drift lines 
depend on the energy level as well. 

(ii) As soon as magnetic lines are defined, the movement along them is defined 
by the quantum number _/' and the energy level r. Moreover, the direction 
of the movement depends on (2.6). In contrast, in the two-dimensional case 
the choice of _/' is automatic and the movement is defined only by the energy 
level which should be fixed. System (1.5) is a dynamical system while (2.1) 
only describes geometrical shape. 

(iii) The movement along magnetic drift hnes is given in scaled time; the 
real time is fit. According to the theory of sections 2.2 and 2.4, t should 
not be too large. Moreover, under condition 

(2.7) |D"(^)|<cz/ \/a:l<\a\<K, 

it is sufficient to assume that ut is not too large (because we can divide 
operator by u and then the scaled time will be ut). 

(iv) Vf~^ is constant along magnetic drift lines and for big enough /x (e.g., 
H > /7^2) one need only consider lines with —V/{fihF) G 2Z+ + 1. 



3 Supersharp estimates as c/ = 3 
3 . 1 Framework 

Now instead of escape conditions we assume that all trajectories of U^^.t are 
trapped either because all magnetic lines are periodic or by the growing 
potential. We discuss only the latter case but the former could be analyzed 
in the same way. Idea is that magnetic lines drift with a speed 0{fi~^) and 
because of this drift we can increase T* from x 1 to x yu. 

So, let us consider coordinate system in which F^ = F^ = 0, F^ j^ and 

V 

(3.1) F ^ ° ^^ ' ' ^ ■^' '^^' - ^° 

where x' = (xi, X2). 
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Then V/F reaches its minimum z°(x') with respect to X3 at [— Q, Co]; 
we assume that this minimum is unique and non-degenerate we also assume 
that it is unique and moreover 

(3.2) eo< {x3-z\x'))-'d,,j<Co. 

We also assume first that |V( \//F)| > e. Then we can extend time 7"^, = h^^^ 
to T* = e but now we will try to go beyond using canonical form. First we 
have a fast magnetron movement in (xs.^a), then an oscillatory movement 
in (x2, .^2) and finally a slow magnetic drift movement in (xi, ,^1). 

Let us denote by T{x,r) the oscillation time with respect to (x2,^2) 
on energy level 0. Then the drift velocity is n~^H'i^ _,_ji^^^ with bo = V o v|/, 
61 = F o vj/ where H means that we consider Hamiltonian field only with 
respect to (xi,^i). 

Therefore (xi, ^i)-shift will be 

(3.3) ,-/ AVU.,*^.-'/_ ^=1%"- 

which in turn is equal to —2fi^^H' + 0{fi^^) with 

(3.4) r]{x,,^,,r)= I {-bo-b,rydx2 

where z^ = z^(xi, ,^1, r) are roots of 

(3.5) bo(xi,ei,z) + bi(xi,ei,z)r = 0, 
z' < z^ < z+, and 

(3.6) We use the same notation for functions of (xi, ,^1) and for functions of 
x' = (xi,X2) as (xi,X2) = M>(xi,^i), 

(3.7) r<r(xi,6) = -^(xi,6,z°(xi,6)) 
and 

(3.8) r(x„ei,r)= r , ^""^ , 

Jz- V-Do - bir 

is a period of oscillations. Note that due to (3.2) 
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(3.9) \Hri\ X 1 as z^*" are close to z°(xi,^i) 

and therefore shift is observable: /j,^-^ ^ fi^^h. The same is true for any r 
provided \Vri\ x 1 for all r < r(xi,{i). 
Note that we can rewrite this condition 

(3.10) |Vx'r7(x , r)\ > cq Vx , r < r(x ) 
with 



(3.11) r/(x',r)= / {-V-Fr)'ds 

Je{x') 

where x' coordinates on the surface transversal to magnetic lines and ds is 
an element of length along such lines where i{x' , r) indicates the segment of 
the magnetic line ^t{x') where V + rF < 0. Our "new" r] is connected with 
"old" f] in the obvious way. 

Therefore 7"* x 1 is upgraded to T* x /x. 



3.2 Case fih < 1 

The fact that the standard under condition (3.15) below T* x 1 can 
be replaced by T* x /i without adding new singularities implies that the 
standard remainder estimate R""" < C/?^^ and R^^ < Ch^^ could be upgraded 
to R""" < Cfi~^h~^ and R^^ < C^^^h~^ respectively as /i/7 < 1. So we arrive 
to 

Theorem 3.1. Let conditions (0.1)-(0.5), (0.6), 

(3.12) |\/| + |V\/| >eo VxG 6(0,1); 
and 

(3.13) |V-|>eo 

he fulfilled. 

Further, let conditions (3.1), (3.2) and (3.10) he fulfilled. Finally let 
non- degeneracy condition 

(3.14) |V(-)|<e =^ |detHess(— )| >e 
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be fulfilled. Then as fih < 1 estimate 

(3.15) R^^ < C^-'h-^ 
holds. 

Proof. We leave to the reader an easy exercise to prove that aU extra terms 
in R^w pQ^i^ ^g dropped so R^^ becomes R^^. D 

This theorem is illustrated by example 3.5(i). On the other hand, assume 
that 1] has a critical point. Let us introduce 7-admissible partition of unity 
with respect to (xi,^i) with 

(3.16) 7 = 7(xi,,^i, r) := max(e|Vr7|, C7). 

Then the shift during period is of magnitude /i~"^7 and uncertainty principle 
requires fj,~^'j x 7 ^ jj,~^h or 7 ^ hi. Therefore selecting 7 = h2~^ we 
conclude that the contribution of the zone {7 > 7} to R^"^ does not exceed 
Cfi^-^h~^. Meanwhile contribution of the zone {7 < 7} does not exceed 
C/)~^ mes({x : |Vr7(x)| < 7}) which is C/?~^7^ under assumption (3.21) 
below. Then we arrive to the remainder estimate 

(3.17) R^^ < Ci2-'h-^ + a'-' 

with an arbitrarily small exponent 6 > 0. 

One can improve it by using logarithmic uncertainty principle; then we 
can select 7 = C{h\ log h\)2 leading us to the remainder estimate 

(3.18) R^^ < Cfi'^h-^ + Ch-^\ log h\ 

However we can do better than this. Let us select 7 = C/72. The following 
problem is easily accessible by our methods: 

Problem 3.2. (i) Prove that as Q' = Q'{xi, ^i-^hDiJ), Q" = Q"{xi, fi~^hDi,J) 
are 7-admissible elements with 7 > C07 then 

(3.19) \F,^,-irXT{t)rQ',,[)jjm < C^h~WT{^y) 

(ii) Then prove that the contribution of 7-element to R^ with T = e does 
not exceed 

(3.20) ch'h'{f^-' + C-y 
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Then we arrive to the following 

Theorem 3.3. Let all conditions of theorem 3.1 except (3.10) he fulfilled. 
Then under condition 

(3.21) |Vx'r7(x', r)| < e =^ | det Hessy r][x' , r)\>e Vx', r < r(x') 

as fih < 1 estimate (3.15) holds. 

Remark 3.4. (i) As cq/j^"^ < /^ < h^^ we need to check conditions (3.10), 

(3.21) only r = [2j + l)/x/7; also condition (3.14) should be checked only as 
\V-{2j + l)^hF\<e- 

(ii) Estimate o{h^^) holds for fih < 1 under weak non-degeneracy condition 

(3.22) mes({x' : \/ri{x', r) = 0}) = as r = (2j + l)fih Vj; 

(iii) Estimate 0(/i~"^/)^^ + /)^2^ holds for jih <1 under intermediate non- 
degeneracy condition 

(3.23) |Vx'??(x', r)\<t ^ I Hess^/ ri{x' , r)\>e 

Vx', r = (2j + l)^/),j G Z+; 

Example 3.5. Consider 

(3.24) a = (6 - ^/iX2)2 + (6 + \lix,f + \/(x). 

(i)As 

(3.25) V = k^xl + lxi-T (r > 0) 

/c > ensures trapping and / ^ ensures condition (3.10). 

(ii) As 

(3.26) V = k^xl + kxl + kxj -T (r > 0) 

/c > ensures trapping and k y^ 0, I2 y^ ensures condition (3.21). 

(iii) In (ii) assume that li = I2 = 1^ > 0; then we have at least two oscillatory 
movements (with respect to (x2,^2) and slow oscillatory drift movement 
with respect to (xi,^i)). Then we can increase T* to T* ^ fi and under 



Chapter 3. Supersharp estimates as d = 3 13 

appropriate non-commensurability condition we can get remainder estimate 
o{i2-^h-^) instead of 0(/i-^/?-2). 

As /i/) — )■ the third periodic movement (with respect to {x^,^^) iii the 
canonical form) plays role and non-commensurability means that at least 
one of two commensurability conditions should be violated. This is fulfilled 
automatically for non-commensurable k and / while for commensurable 
k and / we get one non-commensurability condition to /i. Then we get 
remainder estimate o{fi~^h~^) instead of 0{^~^h~^). 

(iv) Let 

(3.27) V = k^xl-T (r>0). 

Then there is no drift but while fast rotations have frequency 2/i, normal 
oscillation have frequency 2/c and unless they are commensurable we get 
non-periodic movement and as /i/i — )> we get remainder estimate o(/)~^) 
under appropriate assumption to /i as /i/? — )■ +0. 

Problem 3.6. (i) Investigate case when potential traps only trajectories 
with r > 7 and other trajectories escape. One should expect remainder 
estimate R^^ = o{h~^); 

(ii) Investigate the case when all magnetic lines are periodic and V / F is 
constant along them. 

(Hi) Investigate the case when all magnetic lines are periodic and V / F is 
not constant along them (then trajectories of^r.t with r >? are trapped and 
other trajectories are periodic). 

Example 3.7. 

(3.28) a = ei+ il + (6 - /i«(|x |))2 + w{\x'\) - r 

Then F = (-X2, xi, 0)f(|y|)|x'|-i, F = f{\x'\), f{p) = Sptt > and and all 
magnetic lines are circles. 

(i) Then drift speed is dprj[p) with rj = (— i/i/(p) — r\f[p)\) and it is directed 
along X3. Then under appropriate assumptions ensuring (3.10) or (3.14) one 
can derive corresponding asymptotics. 

(ii) As w[p) = —T and f (p) = const there is no drift and again one can derive 
asymptotics o[h^^) under non-commensurability condition as fih — > 0. 
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(a) Oscillation and drift 



(b) Periodic and drift 



Figure 4: Trajectories for S-dimensional model operator: 
a) Operator (3.24)-(3.25) b) Operator (3.28) 



3.3 Schrodinger-Pauli operator as /i/7 > 1 

Consider now Schrodinger-Pauli operator 



(3.29) 



A = J2 Pjg^^Pk +V- l^hF, with Pj = hDj - fiVj 



J.k 



as fih > 1. We assume that modified condition (3.14) 
(3.14)' 



min|— -2j>/?| + |V( — )| < e =^ | det Hess(— ) | > e 



is fulfilled. Then exactly the same arguments as in the proofs of theorems 
3.1, 3.3 lead us to the following 

Theorem 3.8. Consider Schrodinger-Pauli operator (3.29) with fih > 1. 
Let conditions (0.1)-(0.5), (13.2.1), (3.12) and (3.13) be fulfilled. 
Further, let conditions (3.1), (3.2) and (3.14)' be fulfilled. Then 

(i) Under assumption 

(3.10)' |Vx'r7(x', r)| > eo Vx', r < r(x') Vx', r = IJfih, J G Z+ 

estimate 

(3.30) Rj^^ < Ch-^ + C.fih' 
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where R!^^ is due to the perturbation term x h^ in the canonical formP\- 
(ii) Estimate 

(3.31) R^^ < Cfi 
holds under condition 

(3.21)' |Vx'r7(x , r)| < e ^ | det Hess^' v{x' . r)\ > e 

Vx, r = 2Jfih, J G Z+ 

Proof. We leave details to the reader. D 

Remark 3.9. (i) Estimate Ri^^ = o{fih^^) holds under weak non- degeneracy 
condition 

(3.22)' mes({x' : Vri{x' , r) = 0}) = as r = 2i/i/? \/j; 

(ii) Estimate R'^^ = 0(/i/)^2) holds under intermediate non- degeneracy 
condition 

(3.23)' |Vx'r7(x', r)| < e ^ | Hess^/ 77(x', r)| > e 

W,r = 2jfih,j e Z+; 

3.4 Schrodinger-Pauli operator as /i/7 > 1 
reloaded 

While most likely one can improve (3.30) to 

(3.32) Rj^^ < Ch-^ + Cfiexpieh-^) 

we still want to improve it and (3.31) further. To do this we consider 
operator in the reduced form; a,s ^h > C we need to consider only oper- 
ator Ao{xi, HDi;x2, /7D2) with H = ^"^h and we consider it as /i-pseuodo- 
differential operator with an operator-valued symbol a(xi,^i;x2, /JD2) acting 
in the axillary space EI = ^(]Rx2)- Then under condition (3.10) it is micro- 
hyperbolic which leads us to statements (i)-(ii) of the following 

^^ See remark 13.5.13(ii). We should not care about this before as the sharpest 
remainder estimate was 0(/i/)^''"+'^). 
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Theorem 3.10. (i) Let conditions of theorem, 3.8(i) he fulfilled. Then as 
h G (/i^"^, ho) where /?o is small enough constant 

(3.33) R^^ : = 

I /('e(x,x,0)-(27r)"V^"^e(xi,X2;x3,X3;0)V(x)c/x| < Ch'^ 

where e(x'; X3, 3/3; r) is the Schwartz kernel of the spectral projector e(x'; r) of 
operator a(x') = a(xi,X2;x3, /7D3) and as usual we identify objects depending 
on X and on (xi, ,^1, X2) through M/; 

(a) In particular, for ip = ip{x') 



(3.34) R^^ = I ( e(x,x,0)c/x3 - (27r)-V^"'n(^i.^2; 0))^(x)cy^' 




where n(x'; r) = Tr^j e(x'; r) is an eigenvalue counting function of operator 
a(x'). 

Proof. Proof is standard: as the propagation speed with respect to (xi, ,^1) is 
X yU"-*- we can take 7"* x yU and 7"* '^ fi^h with an arbitrarily small exponent 
S > and then we can launch successive approximations with unperturbed 
operator frozen as Xi = yi; then (xi — yi) = 0{fi^^ 7"*). 

We leave details to the reader. Note only that we need h G {n^^, /Jq) 
because actually tj = ri{xi, C,i, h) and to have | V?7| x 1 we need the above 
assumption. D 

Assume now that condition (3.21) is fulfilled instead. 

Theorem 3.11. Let conditions of theorem 3.8(ii) he fulfilled. Let h G 
{fi^^,ho). Then estimate 

(3.35) R^^< C/)-^|log/i| 

with an arbitrarily small exponent holds. 

Proof. Proof is standard. Consider r] = r]{xi,C,i, h) which as /) G {fJ,^^, ho) 
has non-degenerate critical point. 

Then propagation speed with respect to (xi.^i) speed is x /i~"^7 we 
can take T* >^ jj, and T* x | log/i|/?/7^ and then we can launch successive 
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approximations with unperturbed operator frozen as xi = yi; then (xi—yi) = 
0{fi~^'yT^). Here we must assume that 7 > 7 = yU^2 1 log/i| and contribution 
to the remainder of the zone {7 < 7} does not exceed Cfih'-^'j^. 

We leave details to the reader. D 

Example 3.12. (i) Example 3.5(i),(ii) illustrates theorems 3.10, 3.11. 

(ii) In particular, consider example 3.5(ii). Then (as V' = 1) the num- 
ber of eigenvalues below equals to the number of the lattice points 
{(2/ + l)//i-^/), {2j + l)kh} in the triangle {zi > 0, Z2 > 0, z + 1 + Z2 < r}. 
Transition to "Weyl" expression with respect to (xi.^i) means that we 
replace summation with respect to / by integration thus making an error 
0(1) for each j and the total error 0(/j~"^), resulting in the expression 
const ■^/7-iE.(r-(2i + l)/j)^. 

3.5 Introducing boundary 

According to section 13.7 of [Ivr2] as boundary does not break remainder 
estimate 0{h~^ + fih~^) as d = 3. Let us discuss sharper remainder estimates. 

Magnetic billiards must follow magnetic lines; but magnetic lines are not 
billiards and all of them bounce back from the boundary going the same 
path in the opposite direction. So escape condition in the bounded domain 
fails for sure but these magnetic lines drift. 

Consider what happens as Hamiltonian trajectory reflects from the 
boundary. We need to analyze only the following 

Example 3.13. Let 

(3.36) a(x, = ^ [xl + (6 - MXi)2 + ^f) 

and X = {x : X3 > /cxi} i.e. we assume that magnetic lines are transversal 
to ax. 

Then before reflection movement is described by 

(3.37) xi = yU^"^psin(/it) + xi, X2 = jjT^ p cos[jj,t) + X2, X3 = at + X3, 

ii = pcos{iit), ^2 = /i3ci, 6 = 0- 

with constant x, p > 0, a and after reflection movement is described by 
(3.37) with parameters x, p', a' and with t replaced by t' . 
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One can prove easily that 

(3.38) p^ + a^=p'^ + a'^ 
and 

(3.39) x[=xi, ,x!^ = X2-2fi-^lk, a' = a - 21 

with I E [cr — p\k\,a -\- p\k\] the cosine of the incidence angle. 

So, if /c = i.e. F is orthogonal to dX nothing really happens: x' = x, 
p' = p and a' = —a. 

However if /c 7^ i.e. F is not orthogonal to dX, x jumps by —2p^^lk in 
the direction F-*- D TdX and also energy between movement along magnetic 
line and winding redistributes thus magnetic number j which has sense only 
away from the boundary, jumps by 0{{ph)~^k). 

Therefore even calculation of the total jump seems to be a difficult 
problem; most likely there is no consistent movement but just a wobbling. 

Microlocal implications especially as p is close to h~^ seem to be really 
unclear. This leads to an interesting and very challenging 

Problem 3.14. Repeat analysis of subsection 3 and recover remainder esti- 
mate 0{p^^h^^) as ph < 1 under appropriate assumptions: 

(i) Either magnetic drift of trajectories is larger than the possible effect of 
shifting 

(ii) Or movement along magnetic line is bounded by dX only on one side 
and by a growing potential on the other; then jumps are going only in one 
direction and drift direction is disjoint from it. 

(Hi) Consider Schrodinger-P auli operator as ph > 1 and recover remainder 
estimate 0{h^^) under similar assumptions. 

Problem 3.15. In the general case as ph ^ prove remainder estimate 
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